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Abstract 

Starting from a characterization of admissible Cheataev and vako- 
nomic variations in a field theory with constraints we show how the so 
called parametrized variational calculus can help to derive the vakonomic 
and the non-holonomic field equations. We present an example in field 
theory where the non-holonomic method proved to be unphysical. 



1 Introduction 

At least two different procedures to obtain the field equations for a mechanical 
problem with non integrablc constraints on the velocities have been developed. 
They are respectively called the vakonomic and the non-holonomic method and 
are both based on variational principles where a suitable restriction on the set of 
admissible variations is imposed. In the vakonomic (vak) setting the restriction 
arises from geometric considerations, while in the non-holonomic (NH) case it 
is derived from d'Alembert principle. 

The question of which one of the two methods produces equations the solutions 
of which can be physically observed has been extensively studied and it seems 
(see [LM95]) that, at least for a very large class of mechanical constraints, the 
non-holonomic procedure works better. Nevertheless the vakonomic schema 
proved to give interesting results in other frameworks, such as optimal control 
theory (see for example [BBCM03]). 

In field theory, however, the situation is much less clear: both procedures have 
been generalized to provide field equations and Nother currents in some cases 
(see [MPSWOI, FGR04, BFF07] for vak and [BdLdDS02, VCdLMdD05, KV05] 
for NH) , but it is still not evident which one should be better applied in concrete 
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cases. Moreover no fundamental reason justify the NH method since d'Alembert 
principle cannot be formulated. 

Here we aim at contributing to this debate by reformulating both methods 
in terms of parametrized variational calculus: the use of a paramctrization 
sometime helps to hnd field equations without the need of additional variables 
such as Lagrange multipliers. 

We also provide few examples. In particular we find that if we interpret matter 
conservation as a non-integrable constraint in relativistic hydrodynamics, the 
non-holonomic methods give non-physical results (every section satisfying the 
constraint is a solution), while the vakonomic method can be successfully im- 
plemented. In our knowledge this is the first field theory example in which one 
of the two methods has to be rejected, and, surprisingly enough it is exactly the 
one which works in Mechanics. 

2 Constrained field theories 

Let C — ^ M be the configuration bundle whose (global) sections T(C) represent 
the fields (by an abuse of language we will often denote bundles with the same 
label as their total spaces). Let moreover (a; M , y l ) be a fibered coordinate system 
on C. 

Definition 1 A first order Lagrangian on the configuration bundle C AI 
(m = dimM) is a fibered morphism 

L:J X C — ► A m T*M. 

In local coordinates it can be represented as an horizontal to- form L(x tl , y l ,yl l ) ds 
on J C where ds denote the standard local volume form induced on M by the 
coordinates x^. 

Definition 2 Let C — M be the configuration bundle. A constraint of first 
order with codimension r is a submanifold S C J C of codimension r that 
projects onto the whole ofC. 

The constraint be hence expressed by a set of r independent first-order differ- 
ential equations $W(a;' i , y l , yV) = 0. 

Definition 3 A configuration a € L(C) is said to be admissible with respect to 
S if its first jet prolongation lies in S. The space of admissible configurations 
with respect to S is 

L 5 (C) = { ( reL(C)/Im(jV)eS}. 

Definition 4 The set {C, L, S} where C — > M is a configuration bundle, L 
is a Lagrangian on C and S is a constraint, is called a "constrained variational 
problem". 
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2.1 Vak-criticality 



cpt 

Definition 5 Given a compact submanifold D C M, a vak- admissible varia- 
tion (at first order) of an admissible configuration a € Ts(C) is a smooth one 
parameter family of sections {o~t\te\-i ,i[ C r(7r _1 Z?) such that 

1. cr = a\ D 

2. V<G] -1,1[, 0~t\oz) = cr|g£) 

5. In^^jVtl^o) G ^VCflTS. 

In order to check condition 3, one has to verify that the vertical vector field 
V = ■ij 1 (Tt\t=o satisfies the following condition: 



dy i 



-V 1 



dyi 



0. 



(1) 



Definition 6 We say that an admissible section a € Tg(C) is vakonomically 
critical (or vak-critical) for the variational problem {C,L,S} if VD C M for 
any vak- admissible variation {o~t}te]-i,i[ € r(7r _1 .D), we have 



dt I LOjl<Tt 



= 0. 



An equivalent infinitesimal condition is the following 

VD C M, W e VC s. t. Im[/(y o a)} G TS and V\ 9D = 0, 



dy 7 ' 



ya 



dC 



duV a 



ds = 0. 



2.2 Chetaev criticality 
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Definition 7 Given a compact submanifold D C M a Chetaev- admissible vari- 
ation of an admissible configuration a € r,s(C) is a smooth one parameter family 
of sections {crt}te]-i,i[ C r(7r _1 .D) suc/i i/icrf 

7. cr = (t|_d 

«. We]-l,l[, a t \ dD = a\ dD 

3. the vertical vector field j^a t \t=o = V on a with coordinate expression 

a<t)( a ) 

V = V l d l is such that ^—V l = 0. 

dyi 
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Definition 8 We say that an admissible section a € Ts(C) is Chetaev critical 
for the variational problem {C, L, S\ if VD C M for any Chetaev- admissible 
variation {&t}t£]-i,i[ € r(7r _1 L>) we have 



d 

— Loj a t 



= 0. 

t=o 



dt 

An equivalent infinitesimal condition is the following 



c P t d<P (a > 
VD C M, VV e VC s. t. — — V 1 = and V dD = 0, 



dy a dy-^ V 



ds = 0. 



2.3 Integrable constraints 

Definition 9 Let C and B two bundles on the same base M and let {x^,y l ) 
and (a; M , z A ) be two fibered coordinate systems on C and B respectively. Given a 
fibered morphism f : C — ► B projecting onto the identity of M , with coordinate 
expression z A {x^,y l ), its first order jet prolongation j 1 f is a fibered morphism 
j f : J C — > J X B with coordinate expression 

where the operator d^ , called formal derivative, realizes formally the total deriva- 
tive with respect to x^ . 

Theorem 10 Let S £ J X C be a set of integrable constraints linear in the deriva- 
tives, locally expressed as the zero set of the prolongation = daf( a \x, y) = 
of a morphism f : C — > E (E is a vector bundle). With respect to S any 
Chetaev- admissible variation is also vak- admissible and viceversa. 

Proof: The condition for the vertical vector field V = 4iO-t\t=Q to be relative 
to a Chetaev admissible variation is that 

dyl dy l dy l 

On the other hand, for vak-admissibility the following condition is needed 

~W V + ^i dvV =0 d »-& v +s »e? f duV 

and this is equivalent to 

Now obviously (2) implies (3), while if (3) holds then -^jrf^V 1 has to be a 
constant, and being V vanishing at the boundary, (2) holds too. ft 
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Corollary 11 For any constrained variational principle {C,L,S} with inte- 
grable constraint S any vakonomically critical section is Chetaev critical. 



3 Variational calculus with parametrized varia- 
tions 

Definition 12 A parametrization of order 1 and rank 1 of the set of constrained 
variations is a couple (E,¥), where E is a vector bundle E C, while ¥ is a 
fibered morphism (section) 

P : J X C — ► ( J l E)* ®ji c VC. 

If (x^,y a ,e A ) are local fibered coordinates on E and {d a } is the induced fiber- 
wise natural basis of VC, a parameterization of order 1 and rank 1 associates 
to any section y a (x) of C and any section e A (x,y) the section 

[p a A e A +P a A%,e A }d a (4) 

of VC, where p\ and p a A M are functions of {x v , y h {x), d„y b (x)) , while e A de- 
pends on {x v , y b (x)). 

Definition 13 Given a compact submanifold D C M an admissible variation 
of a configuration a G r(C) on D is a smooth one parameter family of sections 
{crt} te i_i j x[ C T(ir~ 1 D) such that 

1. CT = o-\ D 

2. Vf e] - i,i[, at\dD = <y\dD 

3. there exists a section e <E T(E) such that ^jjVtl^o =< Pl^e > oj l o~ 
and {p o a)\dD = 0. 

Definition 14 The set {C,L,P}, where C M is a configuration bundle, 
L is a Lagrangian on C and ¥ is a parameterization of the set of constrained 
variations, is called a "parametrized variational problem". 

Definition 15 We define critical for the parametrized variational problem {C, L, 
those sections of C for which, for any compact D C M and for any admissible 
variation {o~ t } defined on D one has 



d 

jj Loja t 



= 0. 

t=o 



Accordingly, if we use the trivial parametrization P : C — ► VC* ®c VC that 
to any p £ C associates the identity matrix of V P C then the third condition 
becomes empty and we recover free variational calculus. 

For an ordinary variational problem with Lagrangian L = Cds criticality of a 
section of C is equivalent, in local fibered coordinates y a ), to the fact that 
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for any compact D C M and for any V £ T(VC) such that j 1 (V o g)\qd 
one has 

~dC„ it , „ dC 



D 



—V*(x,y(x)) + —d^(x,y(x)) 



dy 



9yt 



ds = 0. 



Explicit calculations (sec [FF03]) show that the above local coordinate expres- 
sions glue together with the neighboring giving rise to the following global one 



VD C M, VV g VC s. t. (V o a)\ dD = 0, 

/ < SL | j x V > oj l a = 
Jd 

where SL is a fibcrcd morphism 

SL : J X C — > {J l VC)* ®ji c A m T*M. 



(5) 



(0) 



To define criticality for first order parametrized variational problems we have 
to restrict variations to those V G T(VC) in (3) with (V o a)\oD = that 
can be obtained through the parametrization from a section e of E satisfying 
j 1 (e o ct)\qd = 0. If (x^^y 11 , e A ) are local fibered coordinates on E and ($>\e a + 
p\ ^dfiE^da is the local representation of < P| .^e > ojV, criticality holds if 
and only if for any compact D C M for any section e with coordinate expression 
e A (x,y) such that both e A (x,y(x)) = and d^s A (x,y(x)) = for all x £ dD, 
we have 



D 



Fir p>r 

^L M e A + P a A + ^-d M s A + P \ «d v s A ) 



ds = 0. (7) 



To set up a characterization of critical sections in terms of a set of differential 
equations let us introduce the following procedure: let us split the integrand of 
(7) into a first summand that factorizes e A (without any derivative) plus the 
total derivative of a second term (a general theorem ensures that this splitting 
is unique; see [BFF07]). To do this, we integrate by parts the derivatives of e 
in the integrand of equation (7). What we get is 



D 



(E A e A + d^V^ + ¥^d u e A )) ds = 



with 



In [BFF07] we have shown that the coefficients E^, and ¥^ are the compo- 
nents of two global morphisms 

E(L,P) : J 3 C — > E* ® c A m T*M 

F(L, P, 7) :J 2 C — ► ( J 1 ^)* ® ji c A m T*M 

and that to the whole procedure can be given a global meaning in terms of 
variational morphisms and global operations between them. The same can also 
be done for higher order Lagrangians and for higher rank and higher order 
parametrizations. 

3.1 Vak-adapted parameterization 

Definition 16 Let C — M be the configuration bundle. A parametrization 

P s : J X C — > (J X EY ®ji c VC 

of the set of constrained variations is said to be vak-adapted to the constraint 
S C J l C if for all e £ r(_E) and a £ Tg(C) the vertical vector field j (< 
P I j e > oj 1 ^) has image in TS. 

To be vak-adapted to a constraint given by the equations $ < -"- ) = one has to 
check that the parametrization with coordinate expression (4) authomatically 
implement condition (1) or, in formula, that we have 

^- [p\e A + V \ »d»e A ] +^d, [p« A e A + p\ % L e A ] = 0. (8) 

Definition 17 A parametrization Pg vak-adapted to a constraint S is said to 
be vak-faithful on a € r<j(C) to S if for all V € VC such that both Im[j 1 (V o 
a)] € TS and (V o o~)\od = hold, there exist a section e € T(E) such that 
< P I j x e > ojV = V" anrf j^e o <r)|a£» = 0. 

The fundamental problem of the existence of a faithful parametrization that is 
vakonomically adapted to a constraint 5* has been studied recently in [GGR06] , 
where a universal faithful parameterization has been found for any constraint 
satisfying certain (quite restrictive) conditions. However we stress that also 
non-faithful parameterizations can be useful for some specific tasks (see Section 
4.3 and Remark 28). 

Proposition 18 Let a G Tg(C) be a vakonomically critical section for the 
constrained variational problem {C, L, S}; then for any adapted parametrization 
P5 the section a is P5 -critical. 
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Proof: We have 



cr G Ts(C) is vak-critical 

t 

* d f 

VD C M, V adm. var. {crj, — L 








VD C M, VV e VC s. t. Hj^V o a)] G TS and cr| aD = 0, 



/ < sl | > ojV = o 



4(a) 

VD C M, Ve G r(S) s. t. j 2 e\ aD = 0, 



/ ( < <S£| j^P' > fe)oj 3 (T = 



E o j 3 cr = 0. 



Step (a) is not an equivalence since there can be admissible infinitesimal varia- 
tions vanishing at the boundary with their derivatives up to the desired order 
that do not come from sections of the bundle of parameters that do vanish on 
the boundary. The last equivalence holds in force of Stokes's theorem, the van- 
ishing of j l+k ~ 1 e on the boundary and the independence of the generators of E. 



Corollary 19 Let er G Ts(C) be an admissible P5 -critical section for the con- 
strained variational problem {C, L, S} and let also Ps be faithful to S on a; then 
a is vakonomically critical for the constrained variational problem {C,L,S}. 

3.2 Chetaev-adapted parameterization 

Definition 20 Let C M be the configuration bundle. A parametrization 



of the set of constrained variations is said to be Chetaev-adapted to the constraint 
S C J X C if for all e G T(E) and a G T S (C) the vertical vector field < P | j Y e > 
oj 1 ^ satisfies condition 3 of Definition 7. 

In coordinates, given the expression (4) for the parametrization the condition 
reads as 



( J 1 ^)* ®jx c VC 



\p a A e A + P a A %e A ]=0 



(9) 



dy\ 
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Definition 21 A parametrization P5 Chetaev- adapted to a constraint S is said 
to be Chetaev-faithful on a G Ts(C) to S if for all V G VC such that both 
d Qyi ) (V o a) 1 = and (V o cr)|a£> = hold there exist a section e G such 
that < P| j 1 ? > oj l o- = V andj l {£o&)\ gD = 0. 

As we did in the previous Section we can prove the following proposition: 

Proposition 22 Let a G rg(C) be a Chetaev critical section for the constrained 
variational problem {C, L, S}. For any Chetaev- adapted parametrization Pg the 
section a is also Ps -critical. 

Corollary 23 Let a G r<j(C) fee an admissible Pg- critical section for the con- 
strained variational problem {C 7 L,S} with respect to a faithful P5; £/ien cr is 
aZso vakonomically critical. 

4 Examples 

In literature very few examples of Lagrangian field theories with constraints are 
present and the question whether the Chetaev or the vakonomic rule produces 
equations whose solutions are physically observed is still open. Let us present 
here two examples: the first is a classic in Mechanics with non-holonomic con- 
straint, while the second, to our knowledge is the first example of a field theory 
where the vakonomic method seems to be preferable to the non-holonomic one 
(the opposite as in Mechanics). 

4.1 A skate on an inclined plane 

This is the model of a skate (or better a knife edge, as called in [BBCM03]) that 
moves on an inclined plane keeping the velocity of its middle point (that is also 
the unique point of contact with the plane, allowing for rotations) parallel to 
the blade (see figure 1). 



Figure 1: A skate on an inclined plane 




The kinematic variables are the coordinates (x, y) of the contact point and the 
direction 9 of the blade (thus the configuration bundle C is R x R 2 x S 1 — ► R). 
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The Lagrangian is 

L = y( i2 + v 2 ) + ~ m 9c«y 

while the constraint S C J 1 C is given by 

isin^ — y cos 9 = 0. 

4.1.1 The non-holonomic setting 

A Chetaev admissible variation for this constraint is a vector field V £ T(VC), 
locally identified by the three components (V X) V y , Vg) with respect to the natural 
base (tjj, J^, that satisfies condition 3 of Definition 7: 

sin 9V X — cos 6*Vy = 0. 

A parametrization of admissible variations can be found by solving the previous 
equation as follows. 

Let us introduce the subbundle K C VC identified as the kernel of the vector 
bundle morphism (V x , V y ,Vs) i — ► V x — V y . It has a two dimensional fiber and 
using in the fiber coordinates (Wi,W2) the fibered immersion in C reads as 
(Wi,W2) 1 — * (Wi,Wi,W2). In formal language a faithful Chctaev-adapted 
parameterization (with zero order and zero rank) is the fibered morphism 

Ps : C — > V*B ® c VC I Vcr £ T(C),VW £ T(VB) 
< P s | (Wi, W 2 ) > o (x(t), y(t),6(t)) = (cos6Wi,sm6W u W2). 

Varying the Lagrangian along this parameterization one get the following first 
variation formula: 

(5L\ j 1 <V S \W >)o J 1 <j = 
= m±(cos6Wi)' — mgeff (sin 6W\ ) + my (sin 9Wi)' + IQWi = 
= -m[x cos 9 + (g cS + y) sin 9] W\ - I9W 2 + m(xcos9 + y sin 9)W\ + I9W 2 

so that the equations of motion are: 

f x cos 9 + (g e s + y) sin 9 = 
1(9 = 



Remark 24 One can check that these equations are exactly the same that can 
be derived from the traditional rule Ej = -\ a ) > where Ej are the Lagrange 
equations of the unconstrained variational problem (see [BBCM03, GR04]. 
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4.1.2 The vakonomic setting 

A vak admissible variation for this constraint is a vector field V € T(VC), locally 
identified by the three components (V x , V y ,Vg) with respect to the natural base 
~5y> m)' sa tisfies (see Definition 16) the following condition: 

sin 8V X — cos 8V y + (cos 8x + sin 8y)Ve = 0. 

A parametrization of admissible variations can be found by solving the previous 
equation. Let us introduce the vector subbundle K C VC identified by Vq = 0. 
In formal language a faithful vak-adapted parameterization (with order and 
rank both equal to 1) if ycosd + xsirtO ^ is the fibered morphism 

P s : 3 X C — ► (JVS)* ®ji c VC/Vae T(C), VV g r(VB) 

<p.u l m.v.) >«»),»(«):«)) - (y-y.. 1 ^;^' ) • 

Varying the Lagrangian along the parameterization gives the following first vari- 
ation formula: 

{8L\j 1 <V S \W > >o j 1 a = 

tA I SU1 8Vx — cos 8V V \ 

= mxV x - mg cf[ Vy + myV y + 10 — , , . 

V y cos 6 + x sin B J 



Integration by parts of the derivatives of the variations leads to the following 
equations of motion: 



mx = I 



y cos y + x sin ( 



my + mg eS = -I — 



6> cos 6> 



(10) 



y cos tf + xsin( 



Remark 25 One can check that these equations are exactly the same that can 
be derived from the Lagrange multiplier traditional rule (see [BBCM03, GR04]/, 
in fact from the variation of the Lagrangian L' = L + \( a }<&( a > one gets the 
following equations 

{x sin — y cos 8 = 
my = (A sin 8)' 
my + mg e ff = — (A cos 8)' 
18 = A (cos 8x + sin 8y) 

and if y cos 8 + x sin 8 ^ one can solve the last one for A and substitute in the 
others to get again equations (10). 
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Remark 26 Let us notice that due to the particularly simple form of the con- 
straint equation we can solve it for 9 in some open subset of the domain getting 
6 = arctan^- . This also is the fundamental reason for which it is so easy in this 
case to find a vak-adapted parameterization. If one now substitutes this expres- 
sion into the Lagrangian, reducing the configuration bundle but increasing the 
order of the Lagrangian, and then one varies it with respect to the independent 
variable [x,y) one gets again equations (10). 

4.1.3 Comparison 

For a comparison of the solution we defer the reader to [BBCM03], where it has 
been shown that for some given trajectories the vakonomic and non-holonomic 
are very much different. For an experimental test of the real observability of 
vak and NH trajectories in a different mechanical system has been carried on in 
[LM95] where the authors found that realistic trajectories obey NH equations. 

4.2 Relativistic hydrodynamics 

Here we present a field theory example where the vakonomic method and the 
non-holonomic one give very different result. In particular we find that the 
non-holonomic theory seems to be non physical (every admissible section is a 
solution). Let us consider a region M of spacetime with metric g a p (here we 
consider it to be fixed, but if we want to study the coupling with gravity the 
formalism can do it as well, see [BFF07]) filled with a barotropic and isocntropic 
fluid. The world lines of the fluid particles and its density describe completely 
the configuration of the system. There are different methods to describe the 
kinematics. We choose to use a vector density J^ds^ physically interpreted as 
follows: let the unit vector field u^(g a p, J M ) = jjr be tangent to the flow lines 

in every point, and the scalar /o(g Q /3, J p ) = -^r= be the density of the fluid. 

The configuration bundle C is thence that of vector densities of weight —1 (the 
transformation rule of J M for a coordinate change x M = x 11 (x M ) with Jacobian 
matrix M is J' 1 = J' 1 det A/ -1 ). The dynamics of the system is ruled by 
the Lagrangian 

L = -^\g~\[p(l + e{p))] ds 

where the scalar e(p) is the internal energy of the fluid from which we can 
derive the pressure P = p 2 ^ 2 . The vector density cannot take arbitrary 
values because of the continuity equation of the fluid 

= 

that needs to hold. It acts as a non integrable constraint S C J l C on the 
derivatives of the fields. 

Some further details on this physical system can be found in [HE73, Sop76, 
KPT79, Tau54]. 
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4.3 The vakonomic setting 

A vak-admissible variation of the field J with respect to the constraint S is 
given by a vector field V € T(VC) represented in local coordinates as ^75777, 
whose first jet is tangent to S, that is to say d^V^ = 0. 
A vak-admissible parameterization is given by the morphism 

P s : J X C — ► (J X TM)* ®,ji c VC 

such that Vcr G r(C) and VX e T(TM) 

< F s I j x X > ojV = £ X J»J?— = {d v J»X v - J v d v X» + J"d v X v )- 9 



Varying the Lagrangian along the given parameterization gives the following 
expression 

< 5L\j l X > o ff = -Jw [ ^L^L(d v Ji 1 X v - J v d v X» + J»d v X v ) ds 

dp p 

where p = p (1 + e(p)) and the following identities hold 

dp { dp \ 

p— = p + P and p d v I — ) = d v P. 
dp \dp) 

Integrating by parts the derivatives of X one finds the following field equations 



(«X + ^)V V P + (p + P) u"V^ = 0. (11) 

Remark 27 in literature one can also find a different description of the sys- 
tem where the fundamental fields are three scalars R a (x t *) physically interpreted 
as the labels identifying the ''abstract fluid particle" passing throw the point 
x^ 6 M. The quantities J M , and p are then defined as suitable functions of 
the fundamental fields and their first derivatives that automatically solve the con- 
straint. This description is dynamically equivalent to our one and it performs, as 
in Remark 26, a reduction of the configuration bundle ( three fundamental fields 
instead of four) increasing the order of the Lagrangian in such a way that the 
variations of the fundamental fields automatically preserve the constraint. We 
defer the reader to the following references [Tau54, Sop76, KPT79] for further 
details. 



Remark 28 The parameterization we used is vak adapted to the constraint 
d^J^ = 0, but not faithful to it in general. Varying J M along P5 is equiva- 
lent to drag the flow lines along a vector field X vanishing on the boundary 
of integration and adjusting p so that the constraint is preserved (see [HE73], 
section 3.3, example 4)- Depending on the solution sometimes it is possible to 
find variations of J M tangent to S that cannot be produced by means of Ps and 
a vector field X vanishing on the boundary of integration. An example is drawn 
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in figure 2: the twisting do not affect J M on the boundary of the cylinder, nev- 
ertheless it is generated by a vector field X which does not vanish on the upper 
part of the boundary itself! In formal terms we have that £xJ^\od = also 
if X\dD 7^ 0. Nevertheless, this is exactly what we want to do from the phys- 
ical viewpoint: when we think to the congruence of curves that represent our 
fluid and we imagine to vary them leaving the boundary fixed, we mean that the 
curves are fixed not only their tangent vectors! To support our choice to vary 
fields along our parametrization we stress that also in the alternative approach 
of "abstract fluid particles" variations leave unchanged the particle identifica- 
tion on the boundary. Solutions of the Euler- Lagrange equations (11) are not 
necessarily vak-critical solutions of the variational problem with the constraint 
S given by d^J^ = 0, but anyway they represent the motions of physical fluids. 




4.3.1 The non-holonomic setting 

A Chetaev-admissible variation of the field J with respect to the constraint 
$ = d^ J^ = is given by a vector field V £ T(VC) represented in local 
coordinates as -Jttt that satisfies condition 3 of Definition 7: 

£j-V» = <=► V = 0. (12) 

For the relativistic fluid, thus, the unique Chetaev-admissible variation is iden- 
tically vanishing. One cannot define any non-trivial variational framework with 
an empty set of variations and insisting on this route one obtain that any section 
a € r(C) is Chetaev-critical and thanks to proposition 22 also Ps-critical for ev- 
ery Chetaev adapted parameterization. This conclusion is clearly non-physical. 
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5 Conclusions 

We have shown how the parametrized variational calculus can contribute to the 
study of non-holonomic and vakonomic field theories. We have defined the no- 
tion of vak-criticality and Chetaev-criticality and compared them with the one of 
criticality with respect to a vak-adapted or a Chetaev adapted parametrization. 
We have also shown examples in Mechanics and Field Theory. In particular we 
think that relativistic hydrodynamics is the first case where it has been shown 
that the vakonomic method is preferable to the non-holonomic one (the opposite 
result with respect to mechanics!). We still cannot guess why this occurs, nor 
whether it is a general rule for field theories; nevertheless it seems to us that it 
is an interesting occurrence and it deserve to be further investigated. 
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